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Abstract 

Action of a hallucinogenic substance, psilocybin, on internal time representation was investigated in two double-blind, placebo-controlled 
studies: Experiment 1 with 12 subjects and graded doses, and Experiment 2 with 9 subjects and a very low dose. The task consisted in repeated 
reproductions of time intervals in the range from 1.5 to 5 s. The effects were assessed by parameter k of the ‘dual klepsydra" model of internal time 
representation, fitted to individual response data and intra-individually normalized with respect to initial values. The estimates k were in the same 
order of magnitude as in earlier studies. In both experiments, k was significantly increased by psilocybin at 90 min from the drug intake, indicating 
a higher loss rate of the internal duration representation. These findings are tentatively linked to qualitative alterations of subjective time in altered 
states of consciousness. 

© 2008 Elsevier Ireland Ltd. All rights reserved. 

Keywords: Altered states of consciousness; Dual klepsydra model; Psilocybin; Subjective duration; Time experience; Time reproduction 


Awareness of temporal order of events and the ability to 
discriminate temporal durations are typical for the normal, 
wakeful conscious state, whereas alterations of subjective time 
experience are often observed in so-called ‘altered states of con- 
sciousness’ (ASC) [11,8,4,18]. Therefore, artificially induced 
ASCs provide good experimental models for testing concepts 
and hypotheses concerning ‘subjective time’ and, particularly, 
models of internal time representation. 

Psilocybin (4-phosphoryloxy-ALV-dimethyltryptamine) 
induces psychosis-like states, most likely by agonist action 
at serotonin 5-HT2A/1A receptors [19]. Alterations of time 
experience in psilocybin-induced states were subjectively 
reported [6], and recently objectively measured [28]. Psilocybin 
significantly shortened subjects’ reproduction of temporal 
intervals longer than 2.5 s, impaired their ability to synchronize 
to inter-beat intervals longer than 2 s, and slowed down their 
preferred tapping rate. Working-memory deficits and subjective 
changes in conscious state, including disturbances in subjective 
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time sense, were also observed. These effects are in line with the 
evidence for psilocybin’s action via serotonergic modulations 
of prefrontal cortex activity [2], and with the important role 
of prefrontal cortex in time perception and timing behavior 
[14,9,29], 

The aforementioned results concerning the effects of psilocy- 
bin on duration reproduction [28] were based on the evaluation 
of reproduction effects undertaken separately for each stimu- 
lus duration. An alternative strategy is to study the subject’s 
response function across the whole range of stimulus durations, 
using a suitable parametric model, and assessing differential 
effects in terms of model parameters rather than response times. 
All subject’s data are thus taken into account at once; degrees of 
freedom are determined solely by the model’s parameterization, 
not by a particular experimental design. 

The following analyses rely on the ‘dual klepsydra’ model 
(DKM) of duration reproduction [20,23]. The DKM is based 
on an accumulation principle, similarly to the standard 
pacemaker-counter model, also called ‘internal clock’ model 
(ICM) [1], but there are important differences. The ICM emu- 
lates a digital electronic stop-watch, whereas the DKM is rather 
reminiscent of early ‘analog computers’: the inflows and states 
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Fig. 1 . The dual klepsydra model (DKM) of time interval reproduction (adapted 
from [23]). Inflow-outflow units (IOU) 1 and 2 are allocated for the represen- 
tation of the 1st and 2nd perceived interval, respectively. Unit 1 is filled from 
t = 0 up to t = s (presented duration) and loses the accumulated state for t > s. 
Unit 2 is filled from t = s + w on, until equality of states is detected at time 
t = s + w + r, where r is the reproduction response. Inflow i (vertical gray 
arrow) into IOUs is constant during the respective intervals; outflow —Ky t is 
proportional to the accumulated state. 

are continuous variables, not pulse series or pulse counts. The 
DKM consists of two inflow-outflow units (IOUs), allocated 
for the representation of duration perceived during the ‘presen- 
tation’ phase and the ‘reproduction’ phase, respectively (Fig. 1). 
An IOUs state y evolves in time t according to the differential 
equation 

dy, 

, = ; - Ky t , yo = o, (1) 

dr 

i.e., the unit acts as a ‘leaky accumulator’ 1 of inflow i. There 
is no built-in metric of accumulated states; the states of the 
1st and 2nd accumulator can be only compared; the equality 
of states, yf ^ = yp\ is subjectively experienced as ‘equality of 
durations.’ The DKM is not a computational device; it is rather 
a primitive comparator of two elapsed time intervals, built upon 
elementary biophysical principles. 2 

The dual klepsydra model yields a so-called ‘klepsydraic 
reproduction function’ (KRF) [23], predicting the reproduc- 
tion response r to depend on the stimulus duration .v and the 
delay w between the presentation and the reproduction phase, 
r = krAy, w ), where 

krfly, w) = a: - 1 ln(l + p(l — e _ * rs )e _KU '). (2) 

Here, at is the ‘loss rate’ coefficient as in Eq. (1), and rj = i\/ ij 
is a ratio between the inflows i into the 1st and 2nd accumula- 
tor. Using the same sensory stimulus to mark both intervals, we 
postulate equal inflow rates, i\ = D. With rj = 1 fixed by postu- 
late, and w being constant by experimental design, the form of 
the response function is determined solely by parameter at. The 


1 Hence the name of the model: klepsydra = ancient water-clock. 

2 In a tentative neurobiological interpretation ([21], Appendix A) the state 
variable is a relative proportion of excited cells in a neuronal assembly, driven by 

an excitatory input (‘inflow’). Although there is no direct observational evidence 
for the DKM, a brain imaging study of temporal reproduction [9] supports the 
two accumulators hypothesis. 


partial functions: 

f K (s) = krf(.y, u>\i] = 1) (3) 

constitute a one-parametric family with parameter at > 0, includ- 
ing the ‘chronometrically correct’ reproduction r = s as a 
limiting case for k —> 0. Given n stimulus-response pairs 
(sj, rj) . =1 n , parameter k can be estimated by minimizing the 
weighted sum of squared errors (WLSQ): 

E(k) = Y, SisjXrj - f K (sj)) 2 . (4) 

j=l...n 

The weighting function g compensates for the increase of 
response variance with increasing stimulus duration 5 [22]; we 
routinely use g(s) = s _1 . 

The functions (3) match experimental data with very good 
accuracy, and fit well with the ‘progressive shortening’ phe- 
nomenon: the responses r are typically shorter than presented 
durations s [30,17], and the ratio r/s decreases with increasing 
stimulus duration s (cf. Fig. 3a). This effect, also dubbed ‘sub- 
jective shortening’ of temporal intervals in memory [27], is seen 
as a negative curvature of reproduction curves [5], and clearly 
revealed by r/,v-versus-.v data plots [12] (cf. Fig. 3b). The ICM 
can account for the progressive shortening phenomenon only 
by increasing (effective) pacemaker frequency with time, which 
itself calls for an explanation. The DKM with constant inflows 
accounts for this phenomenon naturally, by virtue of the ‘loss 
term’ —icy in Eq. (1). The estimates of k are typically in the 
range from 0.01 to 0.03 s _1 [23,24]; the higher the value of at, 
the more pronounced the curvature of the reproduction response 
curve. 

Data from two double-blind, placebo-controlled studies are 
reported. Experiment 1 was part of a previously published study 
[28]. Reproduction data from that study are re-analysed here, 
while data from Experiment 2 are presented here for the first 
time. Psilocybin was administered in gelatine capsules (1 mg 
and 5 mg), prepared at the Pharmacy of the Cantonal Hospi- 
tal of Aarau, Switzerland. 3 Lactose placebo was administered 
in capsules of identical appearance. After oral intake, psilocy- 
bin is rapidly transformed to psilocin, the pharmacologically 
active metabolite. First effects of psilocybin are subjectively 
perceived 20-40 min after oral intake, reach maximum intensity 
at about 60-90 min, last for another 60-120 min, and usually 
vanish within 6-8 h after drug intake [7], 

Placebo (P) and two doses of psilocybin were used in Exper- 
iment 1, medium (M: 115 |xg/kg body weight), and high dose 
(H: 250p.g/kg body weight), with N = 12 subjects participat- 
ing in the study (six men, six women, mean age 26.8 years). 
Subjects were tested on three separate days (for each condi- 
tion) with an interval of at least 14 days in between; P, M and 
H doses were administered in a counter-balanced order. Perfor- 
mance in the timing task was assessed just prior to administration 


3 Administration of psilocybin was authorized by the Swiss Federal Office of 
Public Health, Bern. The experimental procedures were approved by the Ethics 
Committee of the University Hospital Zurich. Subjects had signed an informed 
consent, and were reimbursed for their participation in the study. 
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of drug/placebo (baseline measures), at 90 min during the antic- 
ipated peak effects, and 240 min after drug intake, when effects 
have decreased substantially (for more details see [28]). 

Placebo (P) and a very low dose of psilocybin (VL: 12 jjig/kg 
body weight) were used in Experiment 2, with N = 9 partic- 
ipants (five men, four women, mean age 48.2 years) and one 
measurement at T = 90 min. The two sessions for the respec- 
tive conditions were separated by at least 14 days and included 
the same psychometric tests as employed in Experiment 1 [28]. 
In this low dose study, only few volunteers noticed a difference 
in state of consciousness or mood between placebo and psilo- 
cybin as measured with the Altered States of Consciousness 
Rating Scale and the Adjective Mood Rating Scale (AMRS), 
the only observed effect being an increase in the AMRS scale 
‘introversion.’ 

The duration reproduction task (DRT) followed the same 
design in both experiments. A 500 Hz tone was presented to 
subjects via headphones for a defined duration s (presentation 
phase); after a constant interval of silence, w = 2 s, the same tone 
was presented again and the subjects had to switch off the tone 
by pressing a key when the duration of the second tone was sub- 
jectively equal to that of the first tone (reproduction phase). Tone 
durations s — 1.5, 2, 2.5, 4, 4.5, and 5 s were presented in ran- 
dom order, with eight repetitions for each ^-duration. Response 
times r were measured with 1 ms precision and averaged over the 
eight repetitions. An individual DRT data-set for each subject 
and condition thus consisted of six stimulus-average response 
pairs ( sj , rj) =1 6 . For each data-set, parameter k was estimated 
by the WLSQ method (4) [3]; the estimates are indicated by a 
superscripted caret, £. 

The design of Experiment 1 should result in 3 doses x 
3 time points x 12 subjects = 108 DRT data-sets. Five data- 
sets of the total of 1 08 were unavailable due to a computer failure . 
Nine of the 103 estimates resulted in k — 0, indicating a bad 
KRF-to-data fit. 4 The non-zero estimates were in the range from 
0.0017 to 0.056 s -1 (median « 0.015 s~'), in good agreement 
with estimates from other DRT studies. Because of the failed 
estimates, a total of 14 cells were missing from the data matrix. 
The treatment of the missing values is explained below, together 
with the description of subsequent statistical analyses. Experi- 
ment 2 yielded 2 doses x 1 time point x 9 subjects =18 DRT 
data-sets. All individual data-sets were available, and all esti- 
mates k positive, in the range from 0.0013 to 0.04 s~ 1 (median 
ss 0.02 s _l ), so there were no missing entries; one subject was 
post hoc excluded (see footnote 5). 

To evaluate net time-dependent effects while discarding 
inter-individual differences, each subject’s k s were transformed, 
separately for each drug dose, to intra-individual ratios: 



(5) 


Here kj stands for the estimate k from DRT data measured at 
time T = 0, 90, or 240 min. Therefore, the initial value <70 = 1 


4 In case of a peculiarly shaped individual response curve (positive curvature), 
the error functional (4) attains a minimum at the lowermost boundary of the 
parameter domain, i.e., ic = 0; negative k values are physically meaningless. 


for all subjects in all conditions (not entering the statistics), while 
<790 and <7240 are comparable across subjects. Where one of the 
three values £q, £ 90, £240 was unavailable due to technical failure 
or bad fit (five cases), the missing value was replaced by the 
geometric mean of the two known values. Where two of the three 
values were unavailable (three cases), they were replaced by the 
only one known value; this results in qo = <799 = <7240 = 1, i.e. 
in favor of the null-hypothesis. The same replacement was used 
in one case where all three values were unavailable due to bad 
fit. 

In Experiment 2, DRT measurements were performed only at 
T — 90 min. Therefore, ratios (5) are not applicable, and we can 
only evaluate the drug effects with respect to (w.r.t.) the placebo 
condition, via intra-individual ratios: 


<7vp = 


£ 9 o(VL) 
£90 (P) 


(6) 


The ratios (5), respectively (6), were log-transformed to remove 
the asymmetry of their distribution, and submitted to one- 
sample, two-tailed t - tests against zero (Table 1). The last line 
in Table 1 ( q ) shows anti-log-transformed means = geometric 
mean of individual q’ s; values higher/lower than 1 indicate a 
decrease/increase of £ w.r.t. the initial value £0 or, respectively, 
to the placebo value £90 (P). 

In Experiment 1 at the time of the maximal drug action, 
T — 90 min from intake, the medium dose of psilocybin slightly 
increased the value of £ by ^ 16%, whereas the high dose signif- 
icantly ( P < 0.02) increased the value of £ by ~ 44% (Fig. 2a). 
Due to our conservative treatment of missing values, the mean 
log q was pulled toward 0, so the real effects may be some- 
what larger. No significant effects were found at time T — 240 
min from intake, except for the decrease of k under placebo 
( P ss 0.05). If this is a real effect, it could reflect the progres- 
sive practice in the reproduction task, leading to more precise 
responses and thus counter-acting the ‘subjective shortening’ 
of reproduced times, assessed by parameter k. In Experiment 
2, duration reproduction at time T — 90 min from drug intake 
revealed a significantly higher value of £, by 83%, compared 
to placebo. 5 Interestingly, this effect was observed with a very 
low dose of psilocybin, which barely affected self-reported sub- 
jective experience. 

The effect measures used to evaluate data from Experiments 
1 and 2 have been defined differently and are not directly com- 
parable (see Eqs. (5) and (6)). The two reported studies were 
carried out with groups of participants differing significantly in 
their age (Experiment 1: mean 26.8 years; Experiment 2: mean 
48.2 years). There is evidence for age-dependent decrease of 
5-HT2A receptor expression level [13,16], so that older subjects 
may show higher sensitivity to psilocybin in terms of cognitive 
effects. For all these reasons, we refrain from interpretation of 
different effect sizes between Experiments 1 and 2. 


5 The values shown in Table 1 . column Experiment 2, are based on a subset 
of N = 8 subjects, after removing one subject with q ^ 15.4 as an obvi- 
ous outlier. The entire group’s statistics (N = 9) was: mean log q = 0.365, 
S.D. log# = 0.424, t = +2.581, P = 0.033. The removal of the extreme value 
acts in favor of the null-hypothesis, yet the reported effect remains significant. 
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Table 1 


Summary of psilocybin effects on q ratios in the two reported experiments 


Dose 
T [min] 

Placebo 


Experiment 1 
Medium 


High 


Experiment 2 
Very Low 

90 

90 

240 

90 

240 

90 

240 

Mean log q 

-0.018 

-0.164 

+0.066 

+0.059 

+0.157 

-0.084 

+0.262 

S.D. log q 

0.326 

0.261 

0.194 

0.279 

0.183 

0.259 

0.311 

t 

-0.196 

-2.174 

1.171 

0.733 

+2.971 

-1.123 

+2.383 

P 

ns 

0.052 

ns 

ns 

0.013 

ns 

0.049 

q 

0.958 

0.686 

1.163 

1.146 

1.437 

0.824 

1.828 


The present results clearly indicate the applicability of the 
dual klepsydra model to duration data in the supra-second 
domain. The estimates of the Toss rate’ k are in good agreement 
with our earlier estimates [21,23], in spite of a relatively narrow 
variation range of durations used in the two reported studies. 
Psilocybin, a substance known to strongly influence the subjec- 
tive experience of time, was found to increase the loss rate of 
internal time representation. The effect proved to be dose depen- 
dent in Experiment 1 , where two different doses were compared 
(Fig. 2 b). In Experiment 2, a significant increase under a very 
low dose of psilocybin w.r.t. placebo was found (Fig. 3). The 
DKM parameter k appears to be a fairly sensitive measure of 
effects of psychoactive substances upon the internal represen- 
tation of time. The interpretation of these findings in terms of 
underlying neural mechanisms is not yet clear, as too little is 
known about the sensitivity of the hypothetical ‘accumulators’ 
to environmental or organismic factors. 

The reported findings may help to understand qualitative 
alterations of time experience, such as ‘time standing still’ 
reported from experimentally-induced ASCs [11], mystical or 
ecstatic states (nunc stans), or in psychopathology [15,10]. The 
range of reliable perception of temporal order is determined by 
an interplay between the relaxation time k~ 1 of the DKM accu- 
mulators, and the ‘diffusion time’ given by the signal-to-noise 



Fig. 2. Effects of psilocybin on DKM parameter k in duration reproduction data 
from Experiment 1. Shown are mean log-transformed ratios q = kj/k o (a) for 
placebo (P), medium dose (M, 115 |xg/kg), and high dose (H, 250[xg/kg) of 
psilocybin at times T after drug intake; (b) the same data for the three different 
doses at T = 90 min, demonstrating a dose-dependent effect. Non-linear scale 
of q ratios shown on the right-hand side. 



Presented duration s [sec] 

Fig. 3. Duration reproduction data from Experiment 2. (a) Group averages 
dzlS.D. of individual mean responses r as a function of stimulus duration s. 
Open circles: placebo; filled circles: very low dose of psilocybin (12 |xg/kg). 
Dashed line: ‘chronometrically correct’ response r = s. The curves shown in the 
background are klepsydraic reproduction functions fitted to the data; light gray: 
placebo ( k = 0.0117 s -1 ), dark gray: psilocybin ( k = 0.0212 s -1 ). (b) Group 
average responses relative to the stimulus duration s, showing decrease of the 
r/s ratio with increasing s. Experimental conditions coded as below, section (a). 

ratio of the influx [26]. A simultaneous increase of the influx 
noise and of the loss rate may cause a fusion of the ‘inner’ and 
‘outer’ horizons of duration experience [25], resulting in a total 
loss of cognition of temporal order (“passage of time”). 
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Supplementary data associated with this article can be found, 
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